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摘 要 :研究 了 考虑 宏观 节理 卸载 效应 的 节理 裂隙 岩 体 内 应 力 波 的 能 量 耗 散 规 律 ,综合 分 析 了 细 观 
裂隙 和 具有 不 同 加 印 载 变形 特性 的 宏观 节理 对 应 力 波 传 播 的 影响 。 首 先 , 研 究 了 5 种 波形 入 射 波 
(方形 波 、 半 正弦 波 、 三 角形 波 、 左 三 角形 波 和 右 三 角形 波 ) 在 节理 裂隙 岩 体 内 的 能 量 透射 特性 ;其 
次 ,比较 了 考虑 趣 载 效应 与 不 考虑 孝 载 效应 的 应 力 波 能 量 透射 特性 的 差异 ;最 后 ,讨论 了 入 射 波 波 
) 形 传播 距离 和 入 射 波幅 值 对 能 量 透射 系数 的 影响 ,并 进一步 揭示 了 入 射 波 波形 、 传 播 距 离 和 入 射 
波幅 值 对 考虑 与 不 考虑 外 载 效应 的 能 量 透射 系数 差异 的 影响 。 研 究 结果 表明 :考虑 卸载 效应 的 能 
量 透射 系数 总 是 小 于 不 考虑 卸载 效应 的 能 量 透 射 系 数 。 另 外 ,节理 裂隙 岩 体内 方形 波 的 能 量 迁 射 
系数 最 大 , 右 三 角形 波 的 能 量 透 射 系数 最 小 ;方形 波 的 能 量 透 射 系数 差异 最 小 , 左 三 角形 波 的 能 量 
透射 系数 差异 最 大 。 能 量 透 射 系数 随 着 传播 距离 的 增加 而 减 小 , 随 着 入 射 波 幅 值 的 增加 而 增 大 ;能 
量 透 射 系数 差异 随 着 传播 距离 的 增加 而 增 大 , 随 着 入 射 波幅 值 的 增加 而 减 小 。 

;关键 词 : 印 载 效应 ; 细 观 裂隙; 宏观 节理 ;应 力 波 传播 
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Abstract: The energy dissipation of stress wave propagation through rock masses with microdefects and 
macrojoint considering unloading effects was investigated. The combined effects of microdefects and macro- 
joints with different unloading and loading models on stress wave propagation were considered. First , the 
transmission properties of five waveforms of incident waves (rectangular wave , half-sinusoidal wave , trian- 
gular wave , left-triangular wave and right-triangular wave) propagating through the rock masses with mi- 
crodefects and macrojoint were analyzed. Then ,the differences in the energy transmission properties of 


stress waves propagating through rock masses with microdefects and macrojoint considering unloading 


收 稿 日 期 :2022-08-15 修 回 日 期 :2022-09-22 
基金 项 目 : 国 家 自然 科学 基金 资助 项 目 (No. 12172019 ) ;北京 市 杰出 青年 基金 资助 项 目 ( No. JQ20039) 
通信 作者 : 范 立 峰 。E-mail: fanlifengQ bjut. edu. cn 
引用 格式 : 王 梦 , 范 立 峰 . 秃 载 效应 对 节理 裂隙 岩 体内 应 力 波 能 量 耗 散 影响 研究 [有 站 应 用 力学 学 报 ,2022 ,39 (5) :869-878. 
WANG Meng,FAN Lifeng. The influence of macrojoint unloading effect on the stress waves energy dissipation in rock masses with microdefects 
and macrojoints| J]. Chinese journal of applied mechanics ,2022 ,39 (5 ) 1869-878. 


870 


4 口 工 
村 日 十 
和 四 1 


ChinaXiv& 


应 用 力学 学 报 
effects and those not considering unloading effects were compared. Finally ,the effects of propagation dis- 
tance and incident wave amplitude on the energy transmission coefficient were discussed. The effects of 
propagation distance and incident wave amplitude on the difference of energy transmission coefficient with 
and without considering unloading effect were further revealed. The results show that the energy transmis- 
sion coefficient considering unloading effect is always smaller than the energy transmission coefficient not 
considering unloading effect. In addition, the energy transmission coefficient of the rectangular wave is the 
largest while the energy transmission coefficient of the right-triangular wave is the smallest. The difference 
of energy transmission coefficient of the rectangular wave is the smallest whereas the difference of energy 
transmission coefficient of the left-triangular wave is the largest. The energy transmission coefficient de- 
creases with the increase of propagation distance but increases with the increase of incident wave ampli- 


tude. The difference of energy transmission coefficient increases with the increase of propagation distance 
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and decreases with the increase of incident wave amplitude. 


Key words :unloading effect ; microdefect ; macrojoint;stress wave propagation 


MORBIUBMES T EEG IS ILS BE BU AS VE e A 
面 -例如 细 观 裂 除 和 宏观 节理 ,这 些 不 连续 结构 面 都 
会 溪 响 应 力 波 的 传播 。 细 观 裂 隙 会 导致 应 力 波幅 值 
误 厌 和 波形 弥散 "” ,宏观 节理 会 导致 应 力 波 的 透射 
各 反射 。 由 于 细 观 裂隙 和 宏观 节理 对 应 力 波 传 
授 丽 影响 机 理 不 同 ,难以 建立 统一 的 模型 来 分 析 细 
观 裂隙 和 宏观 节理 对 应 力 波 传播 的 影响 。 此 外 , 宏 
观 着 理 在 动 荷 载 作 用 下 的 加 载 变 形 路 径 和 印 载 变形 
路 径 通 常 是 不 同 的 ,外 载 效应 对 应 力 波 能 量 耗 散 的 
影响 是 不 可 忽略 的 。 研 究 考虑 印 载 效应 的 节理 
裂 耻 岩 体内 应 力 波 的 能 量 耗 散 规律 有 助 于 着 体 工 程 
"Pardi ht ng Ps 

所 目前 , 岩 体 内 应 力 波 的 传播 已 经 得 到 了 广泛 研 
究 忆 通 常 ,位移 不 连续 模型 被 用 于 研究 应 力 波 在 宏 
观 侍 理 处 的 透射 特性 和 反射 特性 " ,该 模型 假设 
宏观 节理 两 侧 的 应 力 是 相等 的 ,而 由 于 宏观 节理 的 
变形 ,宏观 节理 两 侧 的 位 移 是 不 相等 的 5 。 根 据 宏 
观 节理 的 力学 特性 可 将 宏观 节理 的 变形 分 为 线性 变 
形 和 非 线 性 变形 。Schoenberg!5 针对 线性 变形 宏观 
节理 建立 了 线性 位 移 不 连续 模型 ,推导 了 任意 入 射 
角度 下 应 力 波 透射 系数 和 反射 系数 的 解析 解 。 另 
J|, Bandis-Barton ( BB ) 模型 作为 一 个 双 曲 线 函 数 被 
广泛 用 于 描述 宏观 节理 的 非 线 性 变形 057] Zhao 
等 09 基于 BB 模型 建立 了 宏观 节理 两 侧 的 位 移 不 连 
续 关系 ,提出 了 非 线性 位 移 不 连续 模型 ,推导 了 非 线 
性 变形 宏观 节理 处 应 力 波 传播 的 递 推 方程 。 随 后 ， 
Zhao ^5" 研究 了 多 个 平行 非 线性 变形 宏观 节理 对 
应 力 波 传播 的 影响 ,进一步 验证 了 非 线 性 位 移 不 连 
续 模型 的 正确 性 。 天 然 岩 体 中 除了 宏观 节理 , 细 观 


裂隙 也 会 影响 应 力 波 的 传播 “” 。 细 观 裂隙 对 应 力 
波 传播 的 影响 机 理 和 研究 方法 与 宏观 节理 不 同 。 细 
观 裂隙 在 宕 体 中 通常 具有 无 序 和 高 密度 的 特点 。 等 
效 连续 介质 方法 将 细 观 裂隙 岩 体 视 为 连续 均匀 的 介 
质 来 研究 应 力 波 的 传播 ,并 且 为 了 考虑 应 力 波 在 细 
观 裂隙 岩 体 内 传播 的 衰减 特性 和 频率 依赖 特性 , 提 
出 了 等 效 黏 弹性 介质 作为 等 效 连 续 介 质 来 研究 细 观 
裂隙 岩 体 的 动态 力学 特性 。 研 究 等 效 革 弹性 介 
质 中 应 力 波 传播 的 方法 主要 有 特征 线 方法 和 频谱 分 
析 方 法 。 特 征 线 方法 是 一 种 有 限 差 分 方法 ;频谱 分 
析 方 法 是 一 种 傅 里 叶 谐 波 分 析 方 法 。Fan 等 “提出 
了 位 移 不 连续 模型 与 特征 线 相 结 合 的 方法 ,同时 考 
虑 了 细 观 裂 除 和 宏观 节理 对 应 力 波 传播 的 影响 , 研 
究 了 节理 列队 着 体 内 应 力 波 的 传播 特性 。 

在 上 述 研究 中 , 宏观 节理 的 外 载 变形 被 假设 与 
加 载 变形 是 相同 的 。 然 而 ,实验 研究 表明 ,对 于 填充 
了 一 定量 沙子 或 务 士 的 宏观 节理 , RAEE 
加 载 变形 路 径 不 一 致 “ 。 印 载 变形 在 宏观 节理 的 
整个 变形 过 程 中 是 不 可 忽略 的 ,与 加 载 变形 相 比 较 
印 载 变形 表现 出 明显 的 滞后 性 ”。Ma 等 ”提出 
了 一 个 三 相 介 质 模型 来 研究 印 载 行为 对 应 力 波 传播 
的 影响 ,其 中 宏观 节理 的 印 载 变形 被 视 为 是 线性 的 。 
随后 ,引入 了 两 组 BB 模型 分 别 描述 加 载 和 和 缉 载 作 
用 下 宏观 节理 的 变形 特性 ,验证 了 非 线性 加 载 / 
印 载 行为 对 应 力 波 传 播 的 影响 。 人 研究 结果 表明 ,种 
载 效 应 对 岩 体 内 应 力 波 传播 特性 的 影响 是 不 可 忽略 
的 。 然 而 ,由 于 多 种 复杂 因素 的 影响 ,对 考虑 印 载 效 
应 的 节理 裂隙 岩 体内 应 力 波 的 传播 特性 仍 不 其 
了 解 。 
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本 研究 探讨 考虑 印 载 效应 的 节理 裂 际 岩 体 内 应 
力 波 的 能 量 耗 散 规律 。 首 先 ,分 析 了 不 同 波形 的 入 
射 波 在 考虑 邱 载 效应 的 节理 裂 除 岩 体内 的 传播 , 比 
较 了 考虑 印 载 效应 和 不 考虑 伙 载 效应 对 应 力 波 透射 
特性 的 影响 。 然 后 ,讨论 了 入 射 波 波形 传播 距离 和 
入 射 波 幅 值 对 能 量 透 射 系 数 的 影响 ,分 析 了 入 射 波 
波形 传播 距离 和 和 人 射 波幅 值 对 能 量 透 射 系数 差异 
的 影响 。 


1.1 考虑 和 扼 载 效 应 的 分 离 式 三 特征 线 方法 


二 岩 体 内 莹 这 存在 宏观 节理 和 细 观 列队 等 不 连续 
结构 面 , 不 连续 结构 面 影响 岩 体内 应 力 波 的 传播 。 
(a) 为 宏观 节理 对 应 力 波 传播 的 影响 。 其 中 , 宏 
lasst. 图 1(b) 为 细 

隙 对 应 力 波 传播 的 影响 。 由 于 细 观 橡 隙 形状 和 

M eR a oid 
Eni EPA 因此 , 岩 体 中 宏观 
AER EN ARI a. 


透射 波 


透射 波 


ro 


(b) as 
图 1 宏观 节理 和 细 观 裂 际 对 应 力 波 传 播 的 影响 


Fig.1 Effect of macrojoint and microdefect 


裂隙 


on stress wave propagation 
本 研究 同时 考虑 了 岩 体 内 容 观 节理 和 细 观 裂 际 
应 力 波 传 播 的 影响 ,引入 了 宏观 节理 的 加 凶 载 模 
de 
元 。 图 2 为 考虑 外 载 效应 的 节理 裂 队 岩 体内 应 力 波 
传播 的 研究 方案 。 图 2(a) 为 包含 宏观 节理 和 细 观 


王 梦 ， 等 ; 印 载 效应 对 节理 列队 岩 体 内 应 力 波 能 量 耗 散 影响 万 器 3aXiV 
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列队 的 节理 裂隙 岩 体 。 节 理 裂 队 岩 体 被 等 效 为 均匀 
不 连续 模型 如 图 24(b) 所 个。 沟 匀 不 连续 槛 型 中 通 
常 采用 等 效 黏 弹性 模型 研究 细 观 裂 际 岩 体内 应 力 波 
的 传播 特性 ,通过 将 应 力 波 在 等 效 秋 弹性 介质 内 传 
播 的 特性 等 价 为 应 力 波 在 实际 含 细 观 裂隙 的 岩 体 内 
的 传播 特性 ,来 建立 细 观 裂 际 岩 体 与 应 力 波幅 值 衰 
减 和 应 力 波 传 播 速度 的 关系 。 图 2(b) 中 的 等 效 寿 
弹性 模型 由 Maxwell 单元 与 弹簧 并 联 组 成 。 另 外 ， 
采用 位 移 不 连续 模型 研究 宏观 节理 处 应 力 波 的 传播 
特性 。 位 移 不 连续 模型 中 位 移 不 连续 条 件 可 由 宏观 
节理 的 变形 特性 确定 。 在 传统 研究 中 ,宏观 节理 的 
印 载 变形 被 假设 与 加 载 变 形 是 相同 的 , 即 不 考虑 名 
载 效应 对 应 力 波 传播 的 影响 。 在 本 研究 中 ,为 了 考 
虑 仓 载 效应 对 应 力 波 传 播 的 影响 ,引入 了 宏观 节理 
的 加 载 模型 和 仓 载 模型 ,其 中 宏观 节理 的 加 载 变形 
TTE 55 LRL REDE FEES n] , n d 2 (e) Hrs s AA, 
f] 2 (d) P EEE T A IEE RIU 53 AN A TES IBI SALAM 
时 节理 裂 际 岩 体 内 透射 波 的 差异 ,从 图 2(d) 中 可 以 
看 到 和 纯 载 效应 对 透射 波 能 量 的 影响 是 不 可 忽略 的 。 
本 研究 基于 应 力 波 传播 理论 分 析 了 节理 裂隙 岩 体内 
印 载 效应 对 应 力 波 能 量 耗 散 的 影响 。 


考虑 宏观 节理 印 载 变形 
& | | aR $ ^ E 
$ 印 载 模型 /| 效应 SEP 
R © ui 
Al maeu EIE 
Í ; 
H 


宏观 节理 闭合 量 dimm 时 间 Wms 
(c) (d) 
图 2 ZEHRA TERRAE A 
应 力 波 传播 的 研究 方案 


Fig.2 Scheme for the investigation of stress wave 


propagation in rock masses with microdefects and 


macrojoint considering unloading effects 


基于 图 2(b) 中 的 均匀 不 连续 模型 建立 了 包含 
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三 族 特 征 线 的 分 离 式 三 特征 线 方法 5 。 
第 一 族 特征 线 和 对 应 的 特征 相 容 关 系 分 别 为 


dx JE, * Ey uH 
dt E po =Cy (1) 


g -Ee)bEy 
"A PEE EDT li Q 
其 中 :se 为 岩石 质点 应 变 ;c 为 岩石 质点 应 力 ;i 为 应 
力 波 传播 时 间 ;po 为 岩石 密度 jz 为 岩石 质点 速度 jx 
为 沿 着 应 力 波 传 播 路 径 的 笛 卡 尔 坐 标 ;Cv 为 岩石 波 
速 ;By 为 Maxwell 单元 中 弹簧 的 弹性 系数 ; mu 为 
Maxwell 单元 中 黏 壶 的 黏 性 系数 ; 玖 ,为 与 Maxwell 单 
元 并 联 的 弹 得 的 弹性 系数 。 
第 二 族 特征 线 和 对 应 的 特征 相 容 关系 分 别 为 

E, * Eq 


dv 


i dx 

> dr "E (3) 
D. (Gor E,)E, 

dv = - do + d 4 

Mpeg uem 9 


第 三 族 特征 线 和 对 应 的 特征 相 容 关系 分 别 为 
m dx =0 (5) 
o - E,g) Ea 
d EX + e A i (6) 
ON 根据 式 (1) (3) 和 (5) 可 知 ,经 (%,t) 平 面 的 任 
uu 3 条 特征 线 , 即 应 力 波 在 细 观 裂隙 岩 体内 的 
传 攻 可 以 用 具有 三 族 特征 线 的 xr 平面 来 求解 。 另 
外 -为 了 研究 宏观 节理 处 应 力 波 的 传播 ,假设 三 特征 
Hl x = > 为 宏观 节理 的 位 置 , 从 而 建立 分 离 式 三 特 
MAIE, ME 3 所 示 。 从 图 3 中 观察 到 分 离 式 三 
特征 线 方法 由 三 角形 单元 、 萎 形 单元 和 分 离 式 菱形 
单 铝 组成。 根据 已 知 的 边界 条 件 和 初始 条 件 ,可 利 
用 三 角形 单元 来 计算 着 体 左边 界 点 的 应 力 波 ,利用 
萎 形 单元 来 计算 细 观 裂隙 岩 体内 的 应 力 波 ,利用 分 
离 式 菱形 单元 用 来 计算 宏观 节理 处 的 应 力 波 ,通过 
将 宏观 节理 的 加 印 载 模型 引入 到 分 离 式 菱形 单元 中 
来 研究 宏观 节理 的 外 载 效应 对 应 力 波 传播 的 影响 。 

假设 本 研究 中 节理 裂隙 岩 体 的 初始 应 力 、 初 始 
速度 和 初始 应 变 分 别 为 


Y de 


g(x,t) 0, O«xzxL, t=0 (7) 
v(x,1) 20, Oc«xxL, t=0 (8) 
e(x,t) Z0, O«xxL, 1=0 (9) 


其 中 工 为 应 力 波 的 传播 距离 。 
男 外 ,假设 双 尺 度 不 连续 的 左 侧 承受 动 荷 载 , 则 
边界 条 件 为 
o-g'(0,), i=0 (10) 
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1.2. 岩 体 内 左边 界 点 、 细 观 裂隙 岩 体 处 和 宏观 节理 
处 应 力 波 的 传播 


当 应 力 波 在 岩 体 内 传播 时 ,可 以 首先 求解 图 3 
中 OA 特征 线 上 的 应 力 波 传播 , 即 求解 OA 特征 线 上 
的 质点 应 力 、 质 点 速度 和 质点 应 变 ” 。 随 后 ,可 进 
一 步 求解 图 3 中 4A0t 区 域 的 质点 应 力 、 质 点 速度 和 
质点 应 变 。 对 AO 区 域 的 求解 可 分 为 对 左边 界 点 的 
求解 和 对 内 点 的 求解 ,其 中 对 内 点 的 求解 包括 对 细 
观 裂 隙 岩 体 处 和 宏观 节理 处 的 求解 。 

图 3 中 的 左边 界 特征 线 0; 上 的 质点 应 力 可 根 
据 边 界 条 件 式 (10) 得 到 。 男 外 ,左边 界 特征 线 Ot. 上 
的 质点 速度 和 质点 应 变 可 根据 图 3 中 的 三 角形 单元 
求解 。 例 如 左边 界 点 P(0,]) 的 质点 速度 v(0,7) 和 
质点 应 变 e(0,)) 可 根据 三 角形 单元 中 特征 线 PO, 
J) -PCOL, -1) 和 特征 线 P(0,7) -P(0,j -2) 对 应 的 特 
征 相 容 条 件 求 得 。 


宏观 节理 


三 角形 
Hoo, i Z0 P2.) 
ee P RE ! 引入 宏观 
* EXE 
0,1 xod ! 载 模型 
疡 CAL PC72) 


Ns 


图 3 分 离 式 三 特征 线 方法 
Fig.3 Split three characteristic lines method?" 
根据 式 (4) 可 获得 三 角形 单元 中 特征 线 PCO, 
站 )-P(1,7 一 1) 对 应 的 特征 相 容 条 件 为 
v(0,7) -v(1,-1) = 


s Le.) o j-D1- 


[o(1,7-1) -Ee(1, -1)]E 
(E, + Ey)ny 
根据 式 (6) 可 获得 三 角形 单元 中 特征 线 PO, 
D -P(0,j -2) 对 应 的 特征 相 容 条 件 为 
£&(0,) -e(0,j-2) = 


dx — (11) 


1 , 
E «E, L 00: -0(0,7 -2)]+ 


[o(0,7 -2) - E,e(0,j -2)] E, 
(E, * Ey)ny 
由 于 左边 界 特征 线 Ot 上 质点 应 力 oo (0, 门 已 
知 ,将 质点 应 力 o(0,7) 代 入 式 (11) 可 获得 左边 界 点 


2d (12) 
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上 的 质点 速度 "(0,7) ,将 质点 应 力 r(0,J) 代 入 式 
(12) 可 获得 左边 界 点 上 的 质点 应 变 e(0,7) 。 

图 3 中 细 观 裂隙 岩 体 处 的 质点 应 力 、 质 点 速度 
和 质点 应 变 可 根据 葵 形 单元 求解 。 例 如 细 观 裂隙 岩 
体 处 点 P.C) 的 质点 应 力 a (1,4) 、 质 点 速度 v(i,j) 
和 质点 应 变 e (1,7) 可 根据 菱形 单元 中 特征 线 P, (i, 
D-P,G-1,-1) ER Pj(i, 门 -Pj(i+1,j-1) 和 
特征 线 P, i,j) -Pi i,j -2) 对 应 的 特征 相 容 条 件 
求 得 。 

根据 式 (2) 可 获得 菱形 单元 中 特征 线 P, (i,7)- 
Pi(i--1,j--1) 对 应 的 特征 相 容 条 件 为 

v(i) -v(i—-1, -1) = 

LD -eG -Lj- D] s 
s [oG-14-1) Bg TIMES, 
y CE, * Eu) 


(13) 
根据 式 可 获得 次 形 单元 中 特征 线 P i,j) -P,G + 


1G2 1) 对 应 的 特征 相 容 条 件 为 
v7) -v(i +1,j-1) = 


1 2 : 3 
gx -øo(i+1,j-1)] z 


[c(i*1,/-1) - E,e(i £14 - DIE, 
ANI Xx 
bie (E, * Ex) nu 


> (14) 
兴 根 据 式 可 获得 萎 形 单元 中 特征 线 P, (i,j)- 
Qj -2) 对 应 的 特征 相 容 条 件 为 
ej) NL. 
gage -eGi-21- 
[øli j-2) -Ee(i,j -2)] Eg, 
(E, * Eyu)ny 
根据 式 (13 ) 和 式 (14) 可 获得 细 观 裂隙 岩 体 处 
的 质点 应 力 ar (iJ) ,将 质点 应 力 a Gu) LAG) 
或 式 (14) 可 获得 细 观 裂隙 岩 体 处 的 质点 速度 v Ci, 
D ,将 质点 应 力 o(i,j) 代 入 式 (15) 可 获得 细 观 裂隙 
岩 体 处 的 质点 应 变 (i,j) o 
宏观 节理 处 的 应 力 波 传播 ,可 根据 位 移 不 连续 
模型 计算 得 到 。 位 移 不 连续 模型 由 宏观 节理 处 的 应 
力 连 续 条 件 和 位 移 不 连续 条 件 组 成 。 
宏观 节理 处 应 力 连 续 条 件 为 
O (Xi,t)=0 (x1,t)=o(x,t) (16) 
另外 ,根据 宏观 节理 的 加 印 载 模 型 ,可 得 宏观 节 
理 处 的 位 移 不 连续 条 件 为 


"v 


dt (15) 


oem fainaX ive EAF 
reos nite lao IVETTERRTU 73 


u^ (x,t) -u ^ (x,t) 2d" 
Vunl 
T Lol osi / dri 
(17) 
其 中 :uw 为 宏观 节理 处 质点 位 移 ;d 为 宏观 节理 闭合 
量 ;d 为 宏观 节理 初始 闭合 量 ;du 为 宏观 节理 最 大 
允许 闭合 量 ; 态 为 初始 宏观 节理 刚度 ;上 标 ”+ ”表示 
宏观 节理 后 应 力 波 传播 的 参数 ;上 标 “-” 表 示 宏 观 
节理 前 应 力 波 传播 的 参数 ;上 标 “1” 和 “unl” 分 别 表 
示 加 载 模型 和 钊 载 模型 对 应 的 参数 。 
根据 应 力 连 续 条 件 式 (16) .位 移 不 连续 条 件 式 
(17) 以 及 分 离 式 疤 形 单元 中 特征 线 P5 (i,j)-P,(i - 
1 -1) 特征 线 P (i,j)-P (i 1,7] 1) "REGE 
P,* (i j)-P,* (i,j -2) 和 特征 线 P, (1) -Pz (ij - 
2) 对 应 的 特征 相 容 条 件 , 可 求 得 图 3 中 宏观 节理 处 
点 PiJ) 的 质点 应 力 oij) RIEBE vC i,j) RR 
点 应 变 e(i,])。 
根据 式 (2) 可 获得 分 离 式 净 形 单元 中 特征 线 
P; (i) -P(i - 1,j 1) 对 应 的 特征 相 容 条 件 为 
v (i j)-vli-l1,j-1)= 
Ao (ij) -oli -1,j-1)]+ 
Powyv 
[o(i-1,7-1) -Ee(i-1,-1)]Ey 
CE, + Ey) nu 


l/unl 
=d; 


dx 


(18) 
根据 式 (4) I ARIF AN N IN E JE A Te P AF E 2X 
P,* (7)-PG+1J-I) 对 应 的 特征 相 容 条 件 为 
v’ (i) -v(i*1j-1)- 


1 PE i . , 
orna (4) -a(i*1,/-1)] - 
Laer D EECI Ero Hoa 
(E, * Ex) 
(19) 


根据 式 (6) 可 获得 分 离 式 菱形 单元 中 特征 线 
P,* (i,j)-P,* (i,j -2) 对 应 的 特征 相 容 条 件 为 
e` (ij) -e'(ij-2)- 


1 ewa dus ad 
Bam” ()-eo'(ij-2)] + 
*(ij-2)-E,e*(i,j-2)]E 
[oc (i, 7-2) -Ee (i, -2)] vg 


(E, + Eu) 
(20) 
同样 ,根据 式 (6) PI IRIEN EIN E E Hcr p PE 
[EZR P5 (i) -Py 7-2) 对 应 的 特征 相 容 条 件 为 
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e (IJ) -ee (ij-2)- 
1 EUN -—— 
kimi (4g)-e (G-2)]* 
[c (ij-2)-E,e (ij -2)]Ey 
(E, * Eq) zi 
联 立 式 (16) ~ 式 (19) 可 得 宏观 节理 两 侧 的 质 
点 应 力 为 


t (21) 


uC 
ii aT APTE TER 
[eG e Li- D eeG- LIED, 


poCv 
a (i, -2) : ul T" [o |) 
2000) SPP eau 
[eG € 1 71) -Eeli 1 -D]Eg 
E 
(E, + Ey)ny 
WIøe(i-1,j-1)-E,e(i-1,j-1)]Ey | 
> dx} (22 
ES (E, * Ey) i A 


OO 将 式 (22) 代 人 式 (18) ~ (21) 可 分 别 求 得 宏 
Jtr REA E t CREE v^ (i,j) 和 质点 应 变 e (i, 
门 E 以 及 宏观 节理 右 侧 的 质点 速度 v* (i,j) 和 质点 应 
AES (i,j) o £x EL TAKIE HEB BEC ERE 
阶 岩 体内 的 透射 波 。 
二 本 研究 引入 能 量 透 射 系数 也 来 研究 考虑 印 载 效 
应 的 节理 裂隙 岩 体 内 透射 波 的 能 量变 化 ,能 量 透射 
ES MEE] 
I Y is CuriAt 
i ERE EUREN (23) 
之 woCvor4: 
其 中:o* 为 透射 波 ;o 为 人 射 波 ;为 透射 波 的 起 始 
时 疗 ; 刀 为 透射 波 的 终止 时 间 ; 为 人 射 波 的 起 始 时 
间 记 为 人 射 波 的 终止 时 间 。 
另外 ,本 研究 引入 能 量 透射 系数 差异 A7T. 来 研 
究 考虑 仓 载 效应 与 不 考虑 外 载 效应 的 节理 裂隙 岩 体 
内 应 力 波 的 能 量 透射 系数 差异 ,能 量 透 射 系数 差异 
AT 可 表示 为 


万 (24) 


Ep: TAA HE 4 Bou D P ER RE M P I 72 E 
FT GE Ee 3 E CCS TOM AI E E RU WE FO 1 RR R 
岩 体内 应 力 波 的 能 量 透 射 系数 。 


2 研究 结果 


本 节 分 析 了 不 同 波形 入 射 波 垂直 入 射 到 节理 裂 
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际 岩 体 的 能 量 透 射 规律 , 比较 了 考虑 纯 载 效应 与 不 
考虑 印 载 效应 的 节理 裂隙 岩 体内 的 透射 波 。 在 计算 
中 ,等 效 黏 弹性 模型 的 参数 值 分 别 为 E, 227.0 GPa, 
E, =14.0 GPa 和 mw 235.0 MPa . s。 岩 石 的 密度 为 
po =2 500 kg/m 。 加 载 模 型 中 初始 宏观 节理 刚度 和 
最 大 允许 闭合 量 分 别 为 如 =9.901 GPa/m 和 di, = 
0.561 mm。 鲫 载 模 型 中 初始 宏观 节理 刚度 最 大 允许 闭 
合 量 分 别 为 kr 262. 500 GPa/m 和 d?" =0.059 mm, 
男 外 ,以 半 正 弦 波 为 例 说 明 入 射 波 的 振幅 为 Ao = 
0.5 MPa, 入 射 波 的 频率 为 f =200 Hz。 

图 4 为 5 种 波形 入 射 波 在 节理 裂 际 岩 体内 传播 
的 透射 波 。 图 4(a) ~ (e) 中 入 射 波 波形 分 别 为 方形 
波 、 半 正弦 波 \ 三 角形 波 、 左 三 角形 波 和 右 三 角形 波 。 
另外 ,比较 了 考虑 全 载 效应 与 不 考虑 印 载 效 应 的 节 
理 裂 险 岩 体内 的 透射 波 , 并 将 图 4 中 的 和 人 射 波 和 透 
射 波 代入 式 (23 ) 可 分 别 求 得 5 种 波形 入 射 波 在 节理 
裂 际 岩 体内 传播 的 能 量 透 射 系数 。 

从 图 4(a) ~ Ce) 中 可 以 看 出 ,5 种 波形 入 射 波 
在 考虑 缉 载 效应 与 不 考虑 缉 载 效应 的 节理 裂 际 宕 体 
内 传播 ,其 透射 波幅 值 都 会 发 生 衰减 。 对 于 同一 种 
入 射 波 波形 ,考虑 蔓 载 效应 与 不 考虑 撮 载 效应 的 透 
射 波 的 衰减 程度 是 相同 的 ,并 且 透 射 波 波形 是 相似 
的 。 然 而 ,考虑 介 载 效应 的 透射 波 波形 比 不 考虑 御 
载 效 应 的 透射 波 波形 更 罕 。 考 虑 印 载 效应 时 透射 波 
的 能 量 透射 系数 小 于 不 考虑 印 载 效 应 时 透射 波 的 能 
量 透 射 系 数 , 即 考虑 介 载 效应 时 透射 波 的 能 量 损耗 
比 不 考虑 伸 载 效应 时 透射 波 的 能 量 损耗 多 。 

另外 ,从 图 4(a) 和 图 4(d) 可 以 看 到 ,无 论 是 否 
考虑 仓 载 效 应 ,方形 透射 波 和 左 三 角形 透射 波 都 出 
现 了 明显 抛 点 。 不 同 于 方形 透射 波 和 左 三 角形 透射 
波 ,从 图 4(b)、(c) 和 (e) 可 以 看 到 ,无 论 是 否 考 虑 
18] RUM , 半 正 弱 透 射 波 .三 角形 透射 波 和 右 三 角形 
透射 波 都 是 光滑 连续 的 。 其 次 ,通过 比较 方形 波 SE 
EIX 三 角形 波 . 左 三 角形 波 和 右 三 角形 波 的 透射 
波 ,发 现 无 论 是 否 考虑 印 载 效应 ,方形 透射 波 的 能 量 
透射 系数 是 最 大 的 。 说 明了 无 论 是 否 考虑 伸 载 效 
应 ,方形 透射 波 在 传播 过 程 中 的 能 量 消耗 比 其 他 波 
形 的 能 量 消耗 小 。 从 图 4(a) 中 可 以 看 出 ,考虑 印 载 
效应 的 方形 透射 波 的 能 量 透 射 系 数 为 0. 656 3 ,不 考 
虑 色 载 效应 的 方形 透射 波 的 能 量 透 射 系 数 为 0.7302。 
另外 ,无 论 是 否 考虑 仓 载 效应 , 右 三 角形 透射 波 的 能 
量 透 射 系数 是 最 小 的 。 说 明了 无 论 是 否 考 虑 旬 载 效 
应 , 右 三 角形 透射 波 在 传播 过 程 中 的 能 量 消 耗 比 其 
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他 波形 的 能 量 消耗 大 。 从 图 4(e) 中 可 以 看 出 ,考虑 
缉 载 效应 的 右 三 角形 透射 波 的 能 量 透 射 系数 为 
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0. 508 9 ,不 考虑 卸载 效应 的 右 三 角形 透射 波 的 能 量 
透射 系数 为 0.645 6。 
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3c 入 射 波 波形 对 能 量 透射 系数 和 能 量 透 射 系 数 id oc 
Soc 


差异 的 影响 


EMSC 提示 了 人 射流 波形 对 能 能 量 透射 系数 的 
影响 ,比较 了 节理 裂 际 岩 体 中 考虑 印 载 效应 的 能 量 
透射 系数 与 不 考虑 蔓 载 效应 的 能 量 透 射 系数 。 图 5 
(a) 中 入 射流 波形 为 方形 波 、 半 正弦 波 、 三 角形 波 、 
左 三 角形 波 和 右 三 角形 波 。 在 计算 中 ,入 射 波 幅 值 
固定 为 4 20. 5 MPa ,传播 距离 固定 为 了 =5m。 从 
图 5(a) 中 可 以 看 出 ,5 种 波形 入 射 波 中 方形 波 的 能 
量 透 射 系数 最 大 , 右 三 角形 波 的 能 量 透 射 系 数 最 小 。 
另外 ,无论 是 否 考虑 仓 载 效应 ,能量 透 射 系 数 随 方形 
波 \ 半 正弦 波 \ 三 角形 波 \ 左 三 角形 波 和 右 三 角形 波 
是 依次 减 小 的 。 对 于 任意 给 定 的 入 射 波 波形 ,考虑 
pe E 量 透射 系数 总 是 小 于 不 考虑 凶 载 效应 

能 量 透 射 系数 ,说 明了 当 考 虑 纯 载 效应 时 ,应 力 波 
a 6 量 消耗 较 大 。 从 图 5(a) 中 也 可 以 观察 到 , 考 
虑 仓 载 效 应 和 不 考虑 仓 载 效应 的 能 量 透 射 系数 始终 小 
于 1.0, 说 明了 应 力 波 在 传播 过 程 中 能 量 是 衰减 的 。 


0.4 
方形 波 FERE SAWA 左 三 角形 波 右 三 角形 波 


入 射 波 波形 
(a) 入 射 波 波形 对 能 量 透射 系数 的 影 


能 量 透 射 系数 差异 A7./% 


方形 波 HERE SAWE 左 三 角形 波 右 三 角形 波 
入 射 波 波形 
(b) 入 射 波 波形 对 能 量 透射 系数 的 影响 
5 人 射 波 波形 对 透射 波 能 量 的 影响 


Fig.5 Effect of incident waveform on energy 


四 


of transmitted wave 
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图 5(b) 揭 示 了 入 射 波 波形 对 能 量 透 射 系数 差 
异 的 影响 。 能 量 透 射 系数 差异 为 节理 裂隙 宕 体内 考 
虑 缉 载 效应 和 不 考虑 缀 载 效 应 时 能 量 透射 系数 的 差 
异 。 将 图 5(a) 中 考虑 缉 载 效应 的 能 量 透射 系数 和 
不 考虑 邱 载 效应 的 能 量 透 射 系数 代入 式 (24) ,可 求 
得 图 5(b) 中 的 能 量 透射 系数 差异 。 从 图 5(b) 中 可 
以 看 出 ,不 同人 射 波 波形 的 能 量 透射 系数 差异 是 不 
同 的 ,在 5 种 入射 波 波形 中 方形 波 的 能 量 透 射 系 数 
差异 最 小 ,其 能 量 透 射 系数 差异 为 11. 3% , 而 左 三 
角形 波 的 能 量 透射 系数 差异 最 大 ,其 能 量 透 射 系数 
差异 为 28.6% 。 说 明了 和 缉 载 效应 对 方形 波 能 量 透 
射 系数 的 影响 是 最 小 的 ,对 左 三 角形 波 能 量 透 射 系 
数 的 影响 是 最 大 的 。 


32 传播 距离 对 能 量 透 射 系数 和 能 量 透射 系数 差 
异 的 影响 

(ee 

名 四 6(a) 揭示 了 不 同人 射 波幅 值 下 传播 距离 对 


能 晤 透射 系数 的 影响 ,并 比较 了 节理 裂 际 岩 体 中 考 
起 包 载 效应 的 能 量 透射 系数 与 不 考虑 印 载 效应 的 能 
量 透 射 系 数 。 以 半 正 弦 波 为 例 进行 计算 ,和 人 射 波幅 
[ét 分 别 为 0.5 MPa,2.0 MPa 和 3.5 MPa。 从 图 6 
(全 th 可 以 看 出 ,考虑 印 载 效 应 和 不 考虑 印 载 效 应 
的 能 量 透射 系数 随 着 传播 距离 的 增加 而 减 小 ,这 是 
由 等 考虑 到 岩 体 中 细 观 裂隙 对 应 力 波 传播 的 影响 。 
岩 休 中 细 观 裂 阶 密 度 一 定 的 情况 下 ， 等 效 黏 弹性 横 


通 帮 变 ， 此 时 ， 传 播 距离 越 长 ， 岩 体 中 细 观 裂隙 对 
应 盘 波 传播 的 影响 越 大 ， 从 而 导致 能 量 透射 系数 越 
小 经 应 力 波 的 能 量 消耗 越 大 。 另 外 ， 对 于 任意 给 定 
的 四 播 距离， 考虑 印 载 效应 的 能 量 透射 系数 总 是 小 
于 不 考虑 印 载 效应 的 能 量 透 射 系 数 。 从 图 6 (a) 
中 也 可 以 看 出 ， 对 于 任意 给 定 的 传播 距离 ， 无 论 是 
否 考虑 仓 载 效应 ， 和 人 射 波幅 值 越 大 ， 能 量 透 射 系 数 
越 大 。 

图 6(b) 揭 示 了 不 同人 射 波幅 值 下 传播 距离 对 
能 量 透射 系数 差异 的 影响 。 入 射 波幅 值 4 分 别 为 
0. 5 MPa,2.0 MPa 和 3.5 MPa, 将 图 6(a) PH JE HI 
载 效应 的 能 量 透射 系数 和 不 考虑 印 载 效应 的 能 量 透 
射 系数 代入 式 (24) ,可 求 得 图 6(b) 中 的 能 量 透射 系 
数 差异 。 从 图 6(b) 中 可 以 看 出 , 随 着 传播 距离 的 增 
加 ,能 量 透射 系数 差异 线性 增 大 。 另 外 ,对 于 任意 给 
定 的 传播 距离 ,不 同 的 入 射 波幅 值 对 应 不 同 的 能 量 
透射 系数 差异 , 且 入 射 波幅 值 越 大 能 量 透射 系 数 差 
异 越 小 。 


能 量 透射 系数 民 


0 5 10 15 20 25 
传播 距离 L/m 
(a) 传播 距离 对 能 量 透射 系数 的 影响 
一 [一 A=0.5 MPa 


一 一 4-2.0MPa 
28 一 人 一 4=3.53MPa 


0 5 10 15 20 25 
传播 距离 L/m 
(b) 传播 距离 对 能 量 透射 系数 差异 的 影响 
图 6 传播 距离 对 透射 波 能 量 的 影响 


Fig.6 Effect of propagation distance 


on energy of transmitted wave 


3.3 入射 波 幅 值 对 能 量 透 射 系数 和 能 量 透射 系数 
差异 的 影响 


图 7(a) 揭 示 了 不 同 传播 距离 下 入 射 波幅 值 对 
能 量 透射 系数 的 影响 。 其 中 ,传播 距离 工分 别 为 
1.0m,10.0m 和 20.0m。 从 图 7(a) 中 可 以 看 出 , 考 
虑 印 载 效应 和 不 考虑 印 载 效应 的 能 量 透射 系数 随 着 
入 射 波幅 值 的 增加 而 增 大 。 这 是 由 于 非 线 性 变形 宏 
观 节理 的 刚度 随 着 幅 值 的 增加 而 增 大 , 当 入 射 波幅 
值 足 够 大 时 ,宏观 节理 对 应 力 波 传播 的 影响 较 小 ,此 
时 节理 裂隙 岩 体 等 同 于 只 含有 细 观 裂 际 的 宕 体 。 从 
图 7(a) 中 还 可 以 看 出 , 随 着 人 射 波幅 值 的 增加 ,能 
量 透 射 系 数 先 迅速 增加 后 缓慢 增加 。 另 外 ,对 于 任 
意 给 定 的 人 射 波幅 值 ,考虑 伸 载 效应 的 能 量 透 射 系 
数 总 是 小 于 不 考虑 印 载 效应 的 能 量 透 射 系 数 。 对 于 
任意 给 定 的 入射 波幅 值 , 无 论 是 否 考 虑 印 载 效 应 , 传 
播 距离 越 大 ,能 量 透 射 系数 越 小 。 

图 7(b) 揭 示 了 不 同 传播 距离 下 入 射 波幅 值 对 
能 量 透射 系数 差异 的 有 影响。 传播 距离 工分 别 为 
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1.0 m,10.0 m 4120.0 m。 将 图 7(a) 中 考虑 印 载 效 
应 的 能 量 透射 系数 和 不 考虑 仓 载 效应 的 能 量 透 射 系 
数 代 入 式 (24) ,可 求 得 图 7(b) 中 的 能 量 透射 系数 差 
异 。 从 图 7(b) 中 可 以 看 出 ,能 量 透射 系数 差异 随 入 
射 波幅 值 的 增加 而 明显 减 小 。 另 外 ,对 于 任意 给 定 
的 入 射 波幅 值 , 不 同 的 传播 距离 对 应 不 同 的 能 量 透 
射 系数 差异 , 且 传 播 距 离 越 大 能 量 透 射 系 数 差 异 
越 大 。 
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(b) 入 射 波幅 值 对 能 量 透射 系数 差异 的 影响 
图 7 入 射 波幅 值 对 透射 波 能 量 的 影响 


Fig.7 Effect of incident wave amplitude 


chinaXiv 


on energy of transmitted wave 


4 结 论 

本 研究 探讨 了 考虑 印 载 效应 的 节理 裂 际 岩 体内 
应 力 波 传播 的 能 量 耗 散 规律 ,比较 了 不 同人 射 波 波 
形 下 考虑 匠 载 效应 与 不 考虑 伙 载 效应 的 透射 波 ,分 
析 了 入 射 波 波 形 .传播 距离 和 和 人 射 波幅 值 对 能 量 透 
射 系 数 的 影响 ,并 揭示 了 入 射 波 波形 传播 距离 和 和 
射 波 幅 值 对 能 量 透 射 系 数 差异 的 影响 。 得 出 以 下 


结论 。 


1) 将 宏观 节理 的 加 凶 载 模型 引入 到 分 离 式 三 
特征 线 方法 中 ,可 以 有 效 地 研究 逢 载 效 应 对 节理 裂 
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际 岩 体内 应 力 波 传播 的 影响 。 
2) 方形 波 的 能 量 透 射 系 数 最 大 , 右 三 角形 波 的 
能 量 透 射 系数 最 小 。 男 外 ,方形 波 的 能 量 透射 系数 
异 最 小 , 左 三 角形 波 的 能 量 透射 系数 差异 最 大 。 
3) 能 量 透 射 系 数 随 着 传播 距离 的 增加 而 减 小 ， 
随 着 人 射 波幅 值 的 增加 而 增 大 。 另 外 ,能 量 透 射 系 
数 差异 随 着 传播 距离 的 增加 而 增 大 , 随 着 人 射 波幅 
值 的 增加 而 减 小 。 
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